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We derive and investigate an expression for the dynam- 
ically modified decay of states coupled to an arbitrary con- 
tinuum. This expression is universally valid for weak tempo- 
ral perturbations. The resulting insights can serve as useful 
recipes for optimized control of decay and decoherence. 

PACS numbers: 03.65. Ta, 03.65.Xp, 42. 25. Kb, 42.50.Vk 



The quantum Zeno effect (QZE), namely, the inhibi- 
tion of the decay of an unstable state by its (sufficiently 
frequent) projective measurements, has long been con- 
sidered a basic universal feature of quantum systems pi . 
Our general analysis H has revealed the inherent impos- 
sibility of the QZE for a broad class of processes, includ- 
ing spontaneous emission in open space, as opposed to 
the ubiquitous occurrence of the anti-Zeno effect (AZE), 
i.e., decay acceleration by frequent projective measure- 
ments B. Although realistic schemes may well approxi- 
mate such measurements y,^,^ , there is strong incentive 
for raising the question: Are projective measurements the 
most effective way of modifying the decay of an unsta- 
ble state? This question is prompted by two important 
results: (a) A landmark experiment has demonstrated, 
for the first time, both the QZE and AZE by repeated 
on-off switching of the coupling between a nearly bound 
state and the continuum, using cold atoms that are ini- 
tially trapped in an optical-lattice potential 0. (b) It 
has been predicted that periodic coherent pulses, acting 
between the decaying level and an auxiliary one, can ei- 
ther inhibit or accelerate the decay into certain model 
reservoirs M. In both and m, the repeated interrup- 
tion of the "natural" evolution is imperative for decay 
modification. 

In this paper we purport to substantially expand the 
arsenal of decay control, whether measurement-like (i.e., 
accompanied by dephasing) or fully coherent. We de- 
rive a universal form of the decay rate of unstable states 
into any reservoir (continuum) , modified by weak pertur- 
bations with arbitrary time dependence. The results of 
Refs. [||jm,0| are recovered as limiting cases of this uni- 
versal form. Our analysis can serve as a general recipe for 
optimized decay and decoherence suppression for quan- 
tum logic operations ||^ or decay enhancement for the 
control of chaos or chemical reactions [pj . 

Consider the decay of a state |e) via its coupling to a 
system, described by the orthonormal basis {|j)}, which 
forms either a discrete or a continuous spectrum (or a 
mixture thereof). In its most general form, the total 
Hamiltonian is Hq + V{t) + Hi{t), where 



Ho = nuja\e){e\ -t- ?i^ Wj|j)(j|, 



(1) 



with hiOa and titOj being the energies of |e) and |j), re- 
spectively; 



no = E^^^wi^)oi+h.c., 



(2) 



denoting the off-diagonal coupling of |e) with the other 
states, which is dynamically modulated, so as to modify 
the static limit of V effecting the natural decay process; 
and 



H,{t)^hSait)\e){e\+hJ2Sjit)\j){j\, 



(3) 



standing for the adiabatic (diagonal) time-dependent 
perturbations of the energies of the initial (|e)) and fi- 
nal {\j)) states, e.g., AC Stark shifts. 

We write the wave function of the system, with |e) 
populated at i = 0, as 



|*(f)) = a(i)e-*"°*-'/o *"(*')''*' |e) 

-t-E/3,(i)e-'"^*-''i'o^^(*')''*'|j), 



(4) 



the initial condition being |^(0)) = |e). Henceforth we 
treat the generic case wherein the level shifts and the 
temporal modulation of V{t) are independent of j, i.e., 
Sj{t) = Sf{t) and Vje(t) = e{t)fije, e[t) being the modu- 
lation function (Fig. ^ - inset) . Such factorized form of 
the modulation is commonly valid for weak or moderate 
time-dependent fields, which do not appreciably change 
the states of the continuum. One then obtains from the 
Schrodinger equation that the amplitude a{t) obeys the 
exact integro-differential equation |10[| 



^(i)e(i')*(^-^')e"""'*"*''a(t')- (5) 



Here $(t) = fi ^J2j\l^ej\'^e ^^ii*-) is the reservoir re- 
sponse (memory) function and the function e(t) = 

e(i)expH/o<5,/(t')dt'], with 5af{t) = Sa{t) - Sf{t), 
accounts for the modulation of either diagonal or off- 
diagonal elements of the unperturbed Hamiltonian. 

The assumption that the coupling (0) is a weak pertur- 
bation of d3) implies that a{t) varies sufficiently slowly 
with respect to the kernel of Eq. (0), since we then an- 
ticipate [cf. the validity condition (in])] decay rates much 



smaller than the rate of change of the reservoir response 
^{t). One can thus make the approximation a{t') « a{t) 
on the right-hand side (rhs) of Eq. (g). Then one can 
solve Eq. (^ and represent the amplitude modulus of 
level |e) in the form 



\a{t)\ = exp[-i?(i)Q(i)/2], 



(6) 



where we have introduced the fluence Q{t) = /p dT\e{T)\'^, 
and obtained the decay rate in the universal form 

/•oo 

R{t)^2TT dujG{uj+uja)Ft{uj). (7) 



Here G(w) == n-^Re J^ dte"^'<i>{t) = h'^ J^j 



\M 



2j(w- 



ijj) is the coupling spectrum, i.e., the density of states 
weighted by the strength of the coupling to the contin- 
uum or reservoir; Ft{uj) = \et{uj)\'^ /Q{t), with etitL:) = 
(27r)-i/2/Je(i')e"^*'di', is the (normalized to unity) 
spectrum of the modulation function e(i) in the "win- 
dow" (0,t). The result (j^), (0) is valid to all orders of 
t, i.e., it keeps intact the interferences between the mod- 
ulated decay channels and their non-Markovian effects. 
We stress that Eqs. (||), (0) apply to the decay of super- 
posed states J2m'^ni\'^m) (c-g-: in quantum information 
schemes), provided all of them decay and are modulated 
identically. 

We now consider some important consequences of the 
universal form (|g), (0). The modulation spectrum Ft{Ld) 
is roughly characterized by its width vt and the frequency 
shift At = J du;ujFt{uj). A modulation may strongly 
modify the decay rate (analogously to the QZE or AZE) 
whenever Vt + \At\ > £,{uJa), where ^{tOa) is the charac- 
teristic spectral interval over which the weighted density 
of states G{uj) changes near Wq. In particular, if Wa is 
near the edge of the continuum (as for radiative decay 
in photonic crystals or vibrational decay in ion traps, 
molecules and solids), then ^(wa) is the distance between 
uja and the edge j|] (Fig. |l|a). Only in the opposite Hmit, 
i>t + I Ail <^ Ci^a), can one approximately set Ft{Lu) « 
6{uj) in Eq. @, yielding P{t) « exp[-i?GRQ(i)], where 
i?GR = "^T^GluJa) is the extension of the Golden-Rule 
(GR) rate to the case of a time- dependent coupling. 

The modulation function e{t) can be either random 
or regular (coherent) in time, as detailed below. Con- 
sider first the most general coherent amplitude and phase 
modulation (APM) of the quasiperiodic form, e{t) = 
Sfc efce"*"*"*. Here ujk [k = 0, ±1, . . .) are arbitrary dis- 
crete frequencies with the minimum spectral distance 51. 
For a given function e(i) one can obtain —ioJk and e^ as 
the poles and residues, respectively, of the Laplace trans- 
form e{s). If e{t) is periodic with the period H, then 
LDk ~ kfl, and e^ become the Fourier components of e(i). 
For a general quasiperiodic e(t), one obtains 



Q{t)^elt + elY.XkX: 



J.{uJi-uJk)t 



1 



i{uji - uJk) 



(8) 



where e^ = J2k l^fcP equals the average of |e(t)|^ over a 
period of the order of l/fi, Afc = ek/^-c and 

\e,{u;)\'^eltJ2\>^k\'Sivkt/2) 

k 

+4Y.^k\r-^ — ,^_ — '- — . (9) 



k=^l 



27r?7fc?7i 



Here rjk — to — LOk, whereas S{rjkt/2) — 2 sin (rjkt / 2) / irtrij, 
is a sine-function of rjk normalized to 1 . 
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FIG. 1. Decay modification by modulation [Eqs. mfj, 
m)]. Inset: Schematic view of the temporal modulation of 
the shift of level e and its coupling to a continuum, (a) uia 
is near a band edge of G{lj) — Cuj^''^{lj + r)^^6{uj), where 
9{uj) is the unit step function; then a small phase shift (dashed 
peak) is more effective in reducing the decay rate R than large 
phase shifts <j!> ~ vr (dash-dotted peaks) or frequent measure- 
ments/random e(t) (thin curve), (b) Decay rate R (in units of 
_Rgr) in case (a) with uja = CIF: reduction by PM [Eq. (^] 
(curve 1 - (j) = 0.1, curve 2 -</!> = tt) and frequent impulsive 
measurements |2| (curve 3 - QZE) as a function of perturba- 
tion period r (in units of F"^). Curve 1 gives the strongest 
reduction of _R at a given r. 

For t ^ fl^^ the first term on the rhs of {m is a sum of 
peaks, whose spacings are much greater than their width 
2/t. The fast oscillating second term is also peaked at 
Lu = LUk, but we then find that the ratio of the first to 
the second terms, and that of their counterparts in (M), 
is ~ (nt)^^ ^ 1. In the long-time limit, we then ne- 
glect these fast oscillating terms and obtain from Eqs. 
(|)-(|) that P{t) = exp[-R{t)elt], where R{t) in Eq. (|) 
now involves Ft{uj) « J2k \^k\^S{?]kt/'2)- For even longer 
times, exceeding the effective correlation (memory) time 
of the reservoir, tc = maxfc{l/^(ti;Q -\-uJk)}, the functions 
S(r]kt/2) become narrower than the respective charac- 
teristic widths of G(uj) around LUa + cok, and one can set 
S{rikt/2) w S{?]k). Then Eq. (g) is reduced to 



R 



which holds if 



2^^|Afe|2G(L^a+^fe), 
k 



Rtc < 1. 



(10) 



(11) 



This condition is well satisfied in the regime of interest, 
i.e., weak coupling to essentially any reservoir, unless (for 



some k) cua + LOk is extremely close to a sharp feature in 
G{uj), e.g., a band edge [11]. Hence, the long-time limit 
of the general decay rate^) under the APM is a sum of 
the GR rates, corresponding to the resonant frequencies 
shifted by w^, with the weights |Afep. Formula (M) pro- 
vides a simple general recipe for manipulating the decay 
rate by APM. Its powerful generality allows for the opti- 
mized control of decay, not just for a single level, but also 
for hand characterized by a spectral distribution P{uja) 
(e.g., inhomogeneous or vibrational spectrum). We can 
then choose A^ and w^ in Eq. (nO) so as to minimize the 
decay of (||) convoluted with P{uja). The following limits 
of (11^) will be now analyzed. 

(^Monochromatic perturbation: Let e{t) = eoe~"'"-f*. 
Then R = 2TiG{uja + A), where A = 5a f = const is an 
AC Stark shift. In principle, such a shift may drasti- 
cally enhance or suppress R relative to i?GR- It provides 
the maximal variation of R achievable with an exter- 
nal perturbation, since it does not involve any averag- 
ing (smoothing) of G{ijj) incurred by the width of Ft{u)): 
the modified R can even vanish, if the shifted frequency 
Wq -t- A is beyond the cutoff frequency of the coupling, 
where G(w) = (Fig. ||a,b). Conversely, the increase of 
R due to a shift can be much greater than that achiev- 
able with the AZE |g]. In practice, however, AC Stark 
shifts are usually small for (CW) monochromatic per- 
turbations, whence pulsed perturbations should often be 
used, resulting in multiple Wfc shifts according to (|lO| ) . 

(ii) Impulsive phase modulation (PM): Let the phase 
of the coupling amplitude jump by an amount at times 
T, 2t, .... Such modulation can be achieved by a train 
of identical, equidistant, narrow pulses of nonresonant 
radiation, which produce pulsed AC Stark shifts Saf{t) 
in (I). Now e(t) = e*!*/^'^, where [. . .] is the integer 
part. One then obtains that Cc = 1 and Q{t) — t. The 
decay is given by Eqs. (g) and (0), where Ft{uj) can be 
obtained in a closed form. For sufficiently long times 
one can use Eq. (pn). The poles and residues of e(s) == 



(1 

|A 

(f> < 1 



-)/Hi 



^lct>~ 



Afe|2 = 4sin^((/)/2)/(2/c7r 
the k 

|2 



^'^)] yield Wfc — 2kiT /t — (jy/r and 
Y . For small phase shifts, 
peak dominates, |AoP « 1 — 0^/12, 
whereas \\k\^ ~ (fP / Att"^ k'^ for fc 7^ 0. In this case one 
can retain only the fc = term in Eq. (|lC| ) [unless G{u)) 
is changing very fast]. Then the modulation acts as a 
constant shift A = —(P/t. With the increase of |0|, the 
difference between the fc = and k — 1 peak heights 
diminishes, vanishing for ip — ±7r. Then |Aop = |Aip = 
4/7r^, i.e., Ft{uj) for </) = ±7r contains two identical peaks 
symmetrically shifted in opposite directions (Fig. |l|a) [the 
other peaks |Afe|^ decrease with k as (2fc — 1)^^, totaling 
0.19]. 

The above features allow one to adjust the modulation 
parameters for a given scenario to obtain an optimal de- 
crease or increase of R. The PM scheme with a small (p is 
preferable near the continuum edge (Fig. 0a, b), since it 
yields a spectral shift in the required direction (positive 
or negative). The adverse effect of fc 7^ peaks in Ft{Ld) 



then scales as (jP' and hence can be significantly reduced 
by decreasing \(j)\. On the other hand, if oja is near a 
symmetric peak of G{lj), R is reduced more effectively 
for ~ TT, as in Ref. Q], since the main peaks of Ftiuj) 
at loq and loi then shift stronger with t~^ than the peak 
at ojo — ^4>It for <IC 1. 

(iii) Amplitude modulation (AM) oiVcie coupling arises, 
e.g., for radiative-decay modulation due to atomic mo- 
tion through a high-Q cavity or a photonic crystal [Q or 
for atomic tunneling in optical lattices with time- varying 
lattice acceleration [pUl3] . Let the coupling be turned on 
and off periodically, for the time ti and tq — ti , respec- 
tively, i.e., e = 1 for utq < t < utq + ti and e = for 
uTo + Ti <t < {n + 1)to (n = 0, 1, . . .). Now Q{t) in (|) 
is the total time during which the coupling is switched 
on. This case is also covered by Eq. dlG), where the pa- 
rameters are now found to be e^ — ti/tq, ujk — 2fc7r/To, 
|Ao|' = Ti/ro, |Afc|2 = (ri/To)sinc2(fc7rri/ro) (fc ^ 0). 

It is instructive to consider the limit wherein ri <IC tq 
and To is much greater than the correlation time of the 
continuum, i.e., G(w) does not change significantly over 
the spectral intervals (27rfc/ro, 2n(k-\- 1)/to). In this case 
one can approximate the sum (HO) by the integral (m) 
with Ft{uj) « (ri/27r)sinc^(a;Ti/2), characterized by the 
spectral broadening ^ 1/ri (Fig. g - inset). Then Eq. 
(0) for R reduces to that obtained when ideal projective 
measurements are performed at intervals ti [p|. 

Thus the AM scheme can imitate measurement-induced 
(dephasing) effects on quantum dynamics, if the inter- 
ruption intervals tq exceed the correlation time of the 
continuum. This indeed has been observed [0 for atom 
tunneling in optical lattices whose tilt (acceleration) was 
periodically modulated as above. For its analysis we have 
used the approximate expression for <I>(t) obtained in [[l3|, 
which yields the reservoir spectrum G(w -\- uja) (Fig. 2 
- inset), with one maximum at cj ~ cjg, htOg being the 
lattice band gap. The decay probability P{t), calculated 
in Fig. (curves 1-4) for parameters similar to 0] , com- 
pletely coincides with that obtained for ideal impulsive 
measurements at intervals ti [0 and demonstrates either 
the QZE (curve 2) or the AZE (curve 3) behavior. 

The universal Eq. (M), which is a result of unitary 
analysis, is valid also when e{t) is a stationary random 
process. If such a process is characterized by the corre- 
lation time v~^ , one can use a master equation to show 
that, for t ^ i/~^, we have P{t) « e~^*, where the decay 
rate (provided that R <^v) still has the general form (|^), 
but with 



Ft{uj) -^ F{uj) = TT-^e-^Rc/^ e*{t)e{Q)e''^*dt, (12) 

F{ijj) being the normalized spectrum of the random pro- 
cess and e^ = |e(t)|^, where the overbar denotes ensemble 
averaging. Expression (^) with the substitution (|l^) is 
completely analogous to the universal formula describing 
measurement effects on quantum evolution in B]. This 
analogy between unitary and measurement effects stems 
from the ability to emulate projective measurements by 



the dephasing of the level evolution caused by classical 
random fields |§,g. 




FIG. 2. 
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Coupling time ()j.s) 

Tunneling of sodium atoms in optical lattices 
perturbed by AM scheme: the decay probability P{t) as a 
function of the total coupling time. Curves 1,4- decay with- 
out modulation. Curve 2 - QZE (decay slowdown compared 
to curve 1) for n — 0.8 fis, to = 50.8 /is. Curve 3 - AZE (de- 
cay speedup compared to curve 4) for ti — 2 us, to = 52 fis. 
Inset: The coupling spectrum G{uj + LOa) and the scaled mod- 
ulation function T~^F4to{'^) for the conditions of curve 2. 
Here T = ujgd/{na), where a = 15 km/s^ is the lattice accel- 
eration and d = 295 nm is the lattice period, cug = 91 kHz, 
UgT = 2.05 (for curves 1, 2, 5); uj, "" ' " 
(for curves 3, 4). 



g = 116 kHz, ujgT = 3.32 



There may, however, be a notable difference between 
projections and random-field dephasing. Projective mea- 
surements at an effective rate i^, whether impulsive or 
continuous, usually result in a broadened (to a width 
v) modulation function F{uj), without a shift of its cen- 
ter of gravity, A = J dujuiF{uj) « ||,|. This fea- 
ture was shown [p| to be responsible for either the stan- 
dard QZE scaling, R ^ 1/t^, or the AZE scaling. In 
contrast, a weak and broadband chaotic field such that 
\x\I ^ i^B, where / is the mean intensity, i>b is the band- 
width, and X is the effective polarizability, would give rise 
to a Lorentzian dephasing function F(lu) in (12) with a 
substantial shift A = x^- This shift would have much 
stronger effect on R than the QZE or AZE caused by the 

width V ~ x^7 I^B < |A|. 

We have presented here a general theory of dy- 
namically modulated quantum decay, which offers new 
insights into the possibilities of controlling its non- 
Markovian dynamics by off-resonant electromagnetic 
fields. Its unified form (||), (0) encompasses, as special 
cases, all the modulation schemes of current interest, sat- 
isfying the factorization condition [cf. Eq. (^] P,[7|,p^ . 
Whereas its limit (|lj) may imitate measurement effects 



(the QZE and AZE) , the modulation or spectral-shift pa- 
rameters allow us to "engineer" (suppress or enhance) 
more effectively the decay into a given reservoir. Thus, 
measurements are shown to have no advantage as a means 
of either suppressing or enhancing decay compared to 
APM. Moreover, the coherent nature of APM makes it 
much more appropriate than measurements for decoher- 
ence suppression in quantum information applications, 
which require reversible transformations of quantum su- 
perposed states. 
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(A. K.). 
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